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Abstract. We present new empirical and theoretical calibrations of two photometric metal- 
licity indices based on Walraven photometry. The empirical calibration relies on a sample 
of 48 Cepheids for which iron abundances based on high resolution spectra are available in 
the literature. They cover a broad range in metal abundance (-0.5 <[Fe/H]<+0.5) and the 
intrinsic accuracy of the Metallicity Index Color (MIC) relations is better than 0.2 dex. The 
theoretical calibration relies on a homogeneous set of scaled-solar evolutionary tracks for 
intermediate-mass stars and on pulsation predictions concerning the topology of the insta- 
bility strip. The metal content of the adopted evolutionary tracks ranges from Z = 0.001 to 
Z = 0.03 and the intrinsic accuracy of the MIC relations is better than 0.1 dex. 
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1. Introduction 

Classical Cepheids present several advanta- 
geous features when compared with field gi- 
ant stars. These are illustrated by the points 
below: Distance indicators - they obey to 
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well defined optical and near-infrared (NIR) 
Period-Luminosity (PL) relations and their dis- 
tances can be estimated with an accuracy 
of the order of a few percent (Benedict et 
al. 2006; Marconi et al. 2006; Natale et al. 
2007; Fouque et al. 2007). They are also 
the most popular primary distance indicators 
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to calibrate s econdary indicato rs such as the 
SN type la dRiess et al. 2005h . Stellar trac- 
ers - they are distributed across the Galactic 
disk, and therefore, they can be adopted to 
trace the r adial distribution of i ntermediate- 
age stars (lKraft& Schmidt 1963]) . In partic- 
ular, they provide robust constraints on the 
iron and heavy el ements radial gradients 
dLemasle et al. 20071 and references therein). 
Bright stars- they are bright giants that can 
be easily identified due to their intrinsic vari- 
ability and for which accurate photometric and 
spectroscopic data can be collected. However, 
Classical Cepheids also present a few draw- 
backs. Short lifetime - the lifetime they spend 
inside the instability strip is approximately two 
orders of magnitu de shorter than the central 
hydrogen burning dBono et al. 2000) . This to- 
gether with the low spatial density and the 
high reddening account for the limited sam- 
ple of Galactic C epheids curtently known 
dFernie et al. 19951) . Time series - Accurate 
mean magnitudes and colors require a detailed 
time sampling along the pulsation cycle. The 
Cepheid periods range from a few days to hun- 
dred days. Therefore, well sampled multiband 
light curves require long observing runs. In 
this investigation we decide to take advantage 
of the Walraven photometry for 174 Galactic 
Cepheids collected by Pel (1976) and by Lub 
& Pel (1977). The Walraven VBLUW photo- 
metric system was originally designed to de- 
rive the intrinsic properties of early-type stars. 
The bands were selected to measure the fea- 
tures of the hydrogen spectrum, and three of 
the five bands are located in the ultraviolet 
spectral regions (Ay =5405 A, Ab =4280 A, 
Al =3825 A, Au =3630 A, Aw =3240 A). The 
reader interested in a more detailed description 
concerning the ingenious crystal optics filter 
and the instrumentation is refen^ed to Walraven 
& Walraven (1960), Lub & Pel (1977), and 
to Pel & Lub (2007). Because of its sensitiv- 
ity to the Balmer jump and the slopes of the 
B aimer and Paschen continua, the Walraven 
system turned out to be also very useful for 
determining the intrinsic parameters of A, F, 
and G-type stars. For this reason, it has been 
used for t he study o f Galactic an d Magellanic 
Cepheids dPel 1976h : dPel 198l[) and for field 



RR Lyrae stars dLub 19771) . Observations with 
the Walraven five-channel photometer, at- 
tached to the 91 -cm 'Lightcollector' reflec- 
tor, started in 1958 at the Leiden Southern 
Station in Broederstroom, South-Africa. After 
20 years in South- Africa the telescope and the 
photometer were moved to ESO/La Silla in 
Chile. The observations in the new site started 
in March 1979 and continued for a dozen years 
until the decommissioning of the photometer 
in 1991. Thanks to the better photometric con- 
ditions of La Silla, this last phase of the op- 
erational life of the instrument was particu- 
larly fruitful. During 32 years of almost un- 
intertupted observations, a very large amount 
of high-quality data was obtained for many 
types of stars, but the emphasis was on OB- 
associations, clusters. Galactic and Magellanic 
Cepheids (~ 200), Galactic RR Lyrae stars (~ 
100) and faint F/G stars in the disk and in the 
inner Galactic halo. 

2. Data and Theory 

The cmTent Cepheid sar nple was selected frorn 
the observing runs 1962 d Walraven et al. 19641) 
and 1970-1971 ^MMM at the Leiden 
Southern Station (South Africa). The sample 
includes 174 Galactic Cepheids and it is 82% 
complete for all known Cepheids brighter than 
V - 11.0 mag at minimum and south of dec- 
lination -1-15°. Fig. 1 shows the light curves in 
the five bands for three Cepheids of our sam- 
ple, namely FW LUP, RZ VEL and RS PUR 
Data plotted in this figure show that individual 
simultaneous VBLUW measurements present 
an accuracy of the order of few millimag. For 
each object have been collected at least 30 
phase points that properly cover the entire pul- 
sation cycle. This means that the intrinsic ac- 
curacy of the mean magnitudes estimated by 
fitting a cubic spline is of the order of a few 
hundredths of magnitude. 
In order to constrain the possible occurtence of 
systematic errors in the optical photometry we 
are also collecting accurate multiband J, H, K 
NIR data. At present, accurate NIR data are 
available for 65 Cepheids in the Walraven sam- 
ple. For these data the uncertainties of each 
phase point range from 0.005 to 0.007 for 
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Fig. 1. From left to right light curves in 
the Walraven V,B,L,U,W bands for three 
Cepheids. From top to bottom short, interme- 
diate, and long period Cepheids (see labeled 
values). For each band are also plotted the in- 
trinsic scatter of the fit with a cubic spline and 
the amphtude. 



K < 6 mag, deteriorating to about 0.012 at 
K = 8.6. This implies an accuracy in the mean 
magnitudes of about 0.002 - 0.005, depending 
on the number of points. However, the dom- 
inant uncertainty on the mean magnitudes is 
probably due to the absolute calibration which 
is 0.01. For the other Cepheids we use the 
mean NIR magnitudes from the 2MASS cat- 
alog dCutri et al. 2003h . Fig. 2 shows the J, 
H and K-band light curves for four Cepheids 
in our sample. In order to avoid systematic 
uncertainties in the metallicity estimate we 
adopted the reddenings g iven in the Catalog of 
Classical Cepheids dpernie et al. 1995h and re- 
moved the objects that ar e members of binary 
systems (Szabados 2003h . We ended up with a 
sample of 15 1 Cepheids for which we searched 
in the hterature for iron measurements based 
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Fig. 2. From left to right J,H,K-band Ught 
curves for four Galactic Cepheids. The sym- 
bols are the same as in Fig. 1. 



on high resolution spectra. We found accurate 
iron abundances for 48 Cepheids (Andrievsky 
et al. 2002a,b,c; Andrievsky et al. 2004; Luck 
et at. 2006; Mottini et al. 2006; Lemasle et 
at. 2007). Fortunately enough, the calibrating 
Cepheids cover a broad range in metallicity, 
namely -0.5 <[Fe/H]<0.5. 

2.1. Interstellar reddening correction 

Before any physical information can be derived 
from photometric data, all colors and magni- 
tudes have to be corrected for interstellar red- 
dening. In order to constrain the absorption 
coefficients, detailed observations of a sizable 
sample of standard stars in the southern hemi- 
sphere have been performed during three dif- 
ferent periods (1960-61, 1970-71, 1980-81). 
There are very small differences between these 
three absolute calibrations of the Walraven sys- 
tem and transformations are very well stud- 
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ied. The zeropoint was set by the measure- 
ments of one B star passing through the zenith 
in South Africa (w' Sco, spectral type B1.5V, 
V] = 3.96 and (B - V)] = -0.04, HD 144470). 
We underline that we used the data from South 
Africa (1970-71) transformed into the photo- 
metric system as valid for La Silla (1979-1991) 
and then we were able to apply the relations 
between Johnson and Walraven systems. On 
the basis of the quoted data the following re- 
lations were derived by Lub and Pel between 
magnitudes, colors and color excesses in the 
Johnson and in the Wakaven systems: Vj = 
6.886- 2.5V -0.082(V-fi) 
(B-V)j = 2.57(V-B)-1.02(V-B)2 + 0.05(V-B)3 
E{B - V)jIE{V -B) = 2.375 - 0.169(V - B) 
AvjIEiB - V) = 3.17 - 0.16(V - S) - QAlEiy - B) 
while the ratios between different color ex- 
cess in the Walraven system are the following: 
E{B-U)IE{V - B) = OM 
E{U -W)IE{V -B) = QA5 
E(B - L)/E(V -B) = 0.39 
Thanks to the fact that the Wakaven pho- 
tometric system includes five bands, we can 
also define three reddening-free color in- 
dices. For a standard interstellar extinction law 
jCardelli et al. 19891) they are the following: 
[B-U]^ (B - U) - 0.6UV - B) 
[U-W] = (U-W)- 0.45{V - B) 
[B-L] = (B-L)- 0.39(V - B) 
It should be noted that, the Walraven sys- 
tem uses units of logwilntensity) instead of 
magnitudes (-2.5 * logioilntensity)), therefore 
for transforming the distance moduli based on 
Johnson photometry we adopt DMw - -0.4 * 
DMj. 

2.2. Theory 

In order to provide an independent calibra- 
tion of the MIC relation we decided to use 
an homogeneous set of evolutionary mod- 
els for intermediate-mass stars with scaled- 
solar chemical composition. The transforma- 
tion into the observational plane was per- 
formed by adopting a set of atmosphere 
mod els with the same che mical composi- 
tion dCastelhife Kurucz 2003h . The reader in- 
terested in a detailed discussion concerning the 
input physics is referred to Pietrinferni et al. 




Fig. 3. Top: yj,Bj - Vj Johnson Color- 
Magnitude Diagram for calibrating Cepheids 
(diamonds). Solid lines display evolution- 
ary tracks at solar chemical composition 
(Y=0.273, Z=0.0198) and different stellar 
masses (see labeled values). Individual dis- 
tances for calibrating Cepheids were estimated 
using the NIR PL relations by Persson et al. 
(2004) and reddening estimates provided by 
Fernie (1995). The dashed lines show the edges 
of the instability strip. Bottom: Same as the 
top, but for the Walraven y,B-bands. 



(2004). The adopted stellar masses range from 
M = 3.5 - 10.0 Mq, while the metal and he- 
lium content are Z = 0.001, 0.002, 0.004, 
0.008, 0.01, 0.0198, 0.03 and Y = 0.246, 
0.248, 0.251, 0.256, 0.259, 0.273, 0.288. Fig. 
3 shows the comparison between theory and 
observations for both the Johnson Vy, By-bands 
(top) and the Walraven y,B-bands (bottom). 
Individual distances were estimated using NIR 
mean magnitudes and the empirical NIR PL re- 
lations provided by Persson et al. (2004) and 
the reddening corrections provided by Fernie 
(1995). The dashed lines display the predicted 
first overtone blue edge (hotter) and the funda- 
mental red edge (cooler) of t he Cepheid insta- 
bihty strip (Bonoetal. 20051) . 
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Fig. 4. Top: Difference in iron abundance be- 
tween theoretical MIC relations and spectro- 
scopic abundances. The left panel refers to the 
MIC relation based on [L-U] color while the 
right one to the [V-K] color. Bottom: Same as 
the top, but for the empirical MIC relations. 
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Fig. 5. MetalUcity distribution of Galactic 

Cepheids based on theoretical (top panels), 
~ 150 objects, and empirical (bottom pan- 
els), ~ 100 objects. The dashed lines show the 
Gaussian fits, the <t values are also labeled. 



3. Metallicity indices 

In order to provide a new calibration of the 
MIC relation for Cepheids based on Walraven 
photometry we performed a multilinear regres- 
sion fit between the observed (B-L) color, spec- 
troscopic iron abundance and an independent 
color index (CI). In particular, we estimated the 
coefficients of the following MIC relation: 

(B-l) = a+ p[Fe/H] + yCh + S[Fe/H]CIo (1) 

where the symbols have the usual mean- 
ing, and CIq indicates an unreddened color 
index. To estimate the theoretical MIC rela- 
tions we selected, for each set of tracks at 
fixed metallicity (see §2.2), all the points lo- 
cated inside the predicted edges of the in- 
stabihty strip. Then we performed the mul- 
tilinear regression fit (1) between the quoted 



metalhcities and the predicted CIs. We chose 
two different CIs for the multiUnear regression, 
namely (L-U), as originally suggested by Pel & 
Lub (1978), and (V-K) since optical-NIR col- 
ors are good temperature indicators. Moreover 
(L-U) is also a gravity indicator and thus pe- 
riod dependent. By adopting these two col- 
ors we found that for the calibrating Cepheids 
the differences between photometric and spec- 
troscopic abundances were 5[Fe/H](L-ij) = 
0.03, o-(L-u) = 0.19 and 6[Fe/H\v^K) < 
0.01, o-(v-K) = 0.11. For theoretical calibra- 
tion S[Fe/H]^L-u) <0.01, o-(i_c/) = 0.08 and 
6[Fe/H\v-K) =0.01, <T(y-K) = 0.06. In order 
to constrain the possible occurrence of system- 
atic errors introduced by reddening corrections 
we estimated a new set of MIC relations, but 
using reddening free color indices (see §2.1). 
These relations provide the following differ- 
ences: empirical calibrations 6[Fe/H][L-u] = 
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0.02, o-yL^jj} = 0.19 and 5[FelH\v-K] < 
0.01, o-[v-K] - 0.10; theoretical calibrations 
[/.-£/] <0.01, a[L-m = 0.06 and 
6[FelH][v-K] =0.01, a-[v-K\ = 0.06. We un- 
derline that these dispersions are mainly due 
to the standard deviation of the multilinear re- 
gression. In particular, we found crMic(\L-ir\) = 
0.001 and o"m/c([v-^) = 0.006 for theoret- 
ical relations and o-mic{[l-u\) - 0.14 and 
o'Micav-K]) - 0.05 for the empirical MICs. 
Fig. 4 shows the distribution of the difference 
between photometric estimates and spectro- 
scopic iron measurements for theoretical (top 
panels) and empirical (bottom panels) MIC re- 
lations. Note that the left panels display the dif- 
ference for the [L-U] colors while the right one 
for the [V-K] = iV-K)- 2.9(8 - V). 
Finally, we applied the new reddening free em- 
pirical MIC relations to the non calibrating 
Cepheids (~ 100 objects) and estimated their 
individual iron abundances. The same outcome 
applies to the new reddening free theoretical 
MIC relations, but they have been applied to 
the entire sample of ~ 150 Walraven Cepheids. 
Data plotted in Fig. 5 show that the metallicity 
distributions based on theoretical and empiri- 
cal MIC relations based on different color in- 
dices do agree quite well. Moreover and even 
more importantly, the metalUcity range cov- 
ered by the photometric abundances is, within 
empirical and theoretical uncertainties, in good 
agreement with spectroscopic measurements 
(Andrievski et al. 2002; Yong et al. 2006; 
Lemasle et al. 2007) 

4. Discussion 

We provide new empirical and theoretical cali- 
brations of the Walraven metalhcity index. We 
ended up with two independent sets of MIC re- 
lations to estimate Cepheid iron abundances. 
The comparison between iron abundance esti- 
mates based on theoretical and empirical MIC 
relations brings forward two interesting find- 
ings: 

/) Empirical Calibrations - Cepheid iron abun- 
dances based on the empirical calibration agree 
quite well with spectroscopic measurements. 
The difference is on average 0.02 ± 0.02 dex 
with an intrinsic dispersion cr ~ 0.2 dex for the 



relation using the [L-U] reddening free color 
and 0.001+0.002 dex with cr ~ 0.10 dex for the 
relation using the | V-KJ reddening free color. 
ii) Theoretical Calibration - Iron abundances 
based on the theoretical MIC relations agree 
very well with spectroscopic measurements. 
The difference is on average < 0.001 dex with 
an intrinsic dispersion cr = 0.06 dex for the re- 
lation using the [L-U] reddening free color and 
0.01 + 0.02 dex with cr = 0.07 dex for the rela- 
tion using the [V-K] reddening free color. 
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